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A common perception assumes that magnetic 
memories require ferromagnetic materials with a 
non-zero net magnetic moment. However, it has 
been recently proposed that compensated anti- 
ferromagnets with a zero net moment may rep¬ 
resent a viable alternative to ferromagnets.^^® So 
far, experimental research has focused on bistable 
memories in antiferromagnetic metals.^^M^ In the 
present work we demonstrate a multiple-stable 
memory device in epitaxial manganese telluride 
(MnTe) which is an antiferromagnetic counter¬ 
part of common II-VI semiconductors.^^ Favor¬ 
able micromagnetic characteristics of MnTe allow 
us to demonstrate a smoothly varying antiferro¬ 
magnetic anisotropic magnetoresistance (AMR) 
with a harmonic angular dependence on the ap¬ 
plied magnetic field, analogous to ferromagnets. 
The continuously varying AMR provides means 
for the electrical read-out of multiple-stable an¬ 
tiferromagnetic memory states which we set by 
heat-assisted magneto-recording and by chang¬ 
ing the angle of the writing field. We explore 
the dependence of the magnitude of the zero- 
field read-out signal on the strength of the writ¬ 
ing field and demonstrate the robustness of the 
antiferromagnetic memory states against strong 
magnetic field perturbations. We ascribe the 
multiple-stability in our antiferromagnetic mem¬ 
ory to different distributions of domains with the 
Neel vector aligned along one of the three c-plane 
magnetic easy axes in the hexagonal MnTe film. 
The domain redistribution is controlled during 
the heat-assisted recording by the strength and 
angle of the writing field and freezes when suffi¬ 
ciently below the Neel temperature. 

Switching between two memory states in metallic 
antiferromagnets has been demonstrated in FeRh by 
taking advantage of its high-temperature ferromagnetic 
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phase,!^. Similar effects were shown in thin-film IrMn 
tunneling devices utilizing once the exchange-spring with 
an adjacent ferromagnetic layei!^ or by field-coolingP In 
GuMnAs current induced torques were used to switch the 
antiferromagnet even without an auxiliary ferromagnetic 
element and magnetic fields P The two distinct stable 
states of these antiferromagnetic-metal memories showed 
distinct resistances which allowed for their electrical read¬ 
out. The physical mechanism of the read-out is ascribed 
in these devices to the antiferromagnetic AMR. The in¬ 
terpretation appeals to the general principle formulated 
by Louis Neel in his Nobel lecture that: ’’Effects in anti¬ 
ferromagnets depending on the square of the spontaneous 
magnetization should show the same variation as in fer¬ 
romagnetic substances” ^ 

In ferromagnets, a common characteristic feature of 
the non-crystalline AMR is its harmonic cos2(^^ and 
sm2ipB dependence on the angle between the current 
and the applied saturating magnetic field when measured 
with the longitudinal and transverse voltage probes, re¬ 
spectively. In a conventional Hall-bar geometry, the ratio 
of these longitudinal and transverse AMR signals scales 
with the Hall-bar aspect ratio.^ In antiferromagnets, 
however, controlling the Neel order by external magnetic 
field is significantly more difficult which has hindered the 
detection of the characteristic angular dependence of the 
AMR in antiferromagnets. 

Below we present measurements of the harmonic- 
function AMR in our Hall-bar devices made in antifer¬ 
romagnetic semiconductor MnTe and by this directly il¬ 
lustrate the applicability of the above Neel’s principle 
on this basic spintronic phenomenon. Moreover, unlike 
the previously studied bistable antiferromagnetic-metal 
memories, our MnTe devices allow us to set by heat- 
assisted magneto-recording multiple-stable antiferromag¬ 
netic memory states. When removing the writing mag¬ 
netic fields sufficiently below the Neel temperature, the 
harmonic AMR-like signal is still preserved in the zero- 
field read-out traces. As expected for antiferromagnets, 
these multiple-stable states are not erased by even strong 
magnetic-field perturbations as we also demonstrate be- 
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low. 

We use MnTe thin films grown epitaxially on 
InP(lll)A substrates by molecular beam epitaxy. The 
quality of our epilayers is evidenced in Fig. which 
shows a high resolution high angle annular dark field 
(HAADF) image, showing the individual atomic Te and 
Mn columns as bright and dimmer spots , respe ctively. In 
contrast to most thin film MnTe studies^^Hi^, we obtain 
the hexagonal NiAs bulk structur^^ (a-MnTe) due to the 
good lattice matching. X-ray diffraction measurements 
shows that this is the only phase present in our films and 
that c-planes are aligned parallel to the substrate sur¬ 
face. From reciprocal space maps for the 50 nm thick 
epilayer used in the transport experiments we conclude 
that the film is relaxed forming a mosaic structure with 
a block size of 25 ± 5 nm. The relative misorientation 
of the blocks is small with a Gaussian distribution with 
standard deviation of 0.25 ± 0.1°. (For more details on 
the X-ray characterization see Supplementary informa¬ 
tion (SI) Sec. L) 

The semiconducting electronic structure of our a- 
MnTe films is confirmed by optical absorption measure¬ 
ments (see SI Sec. II) in the visible and mid-infrared 
range on epilayers deposited on a transparent SrF 2 (lll) 
substrate. We performed measurements for layers of dif¬ 
ferent thicknesses, which all show the same onset of the 
absorption shown in Fig. ^p. From this we infer an indi¬ 
rect band gap of 1.46 ± Q.I Q eV in agreement with previ¬ 
ously reported bulk values 

Below the Neel temperature, which in bulk a-MnTe 
is 310 kP the magnetic structure consists of ferromag- 
netically ordered Mn planes which are antiferromagneti- 
cally stacke d alon g the c-direction.E^ Neutron diffraction 
experimentand susceptibility measurement (see 
also Fig. Sib) found that the moments lie within the 
Mn-planes, as indicated in Fig. Et- The moderate Neel 
temperature and the expected weak in-plane magneto- 
crystalline anisotropy in the hexagonal structure suggest 
that significant magnetic moment reconfigurations might 
be achievable at moderate applied magnetic fields. Low 
temperature magnetization data measured by the super¬ 
conducting quantum interference device (SQUID) and 
shown in Fig. confirm this expectations. Consistent 
with the antiferromagnetic order, we observe a zero rema¬ 
nent moment and even at high fields the contribution per 
Mn ato m is on ly a minor fraction of it magnetic moment 
of The moderate exchange energy (Neel tem¬ 

perature) and expected weak in-plane anisotropies are 
reflected in the onset of a sizable magnetic moment at 
a moderate threshold of ^ 2 T in the in-plane magnetic 
field sweeps. 

As shown in SI Sec. Ill, the threshold which may be as¬ 
sociated with the spin-flop field in the domains moves to 
lower fields when approaching the Neel temperature. Due 
to the hexagonal crystal structure and antiferromagnetic 
moments oriented in the c-plane, the system can break 
into three types of domains in which the Neel vector 
points along one of the three high-symmetry directions.!^ 


In Of-MnTe single crystals, it was demonstrated that the 
distribution of these domains can be rearranged by cool¬ 
ing the system through the antiferromagnetic transition 
in an applied magnetic field.!^ All these favorable mag¬ 
netic characteristics are the basis of our experiments in 
the MnTe memory devices presented below. 

Transport measurements shown in Fig. [p complete 
the basic characterization of our thin-film a-MnTe sam¬ 
ples. The temperature dependent resistivity shows a 
peak near the Neel temperature associated with crit¬ 
ical scattering off spin fluctuations.!^ Note that criti¬ 
cal anomalies in our films are consistently observed also 
in the susceptibilit 5 ^^ and lattice parameter measure¬ 
ments (see SI Fig. SI, S2). From the Hall effect mea¬ 
surements we obtained low-temperature hole density of 
p = Q X 10^^ cm“^ due to unintentional doping in our 
film and corresponding hole mobility of /i = 43 cm^/Vs. 

We now proceed to the discussion of the antiferromag¬ 
netic AMR and memory functionalities in our MnTe de¬ 
vices. In Figs. [^,b we plot the transverse and longitu¬ 
dinal AMRs, defined as AMR±{ipB) = ^ 

and AMRii{(Pb) = were Rxy{<Pb) and 

^xx(^b) are the transverse and longitudinal resistances 
indicated in the inset of Fig. Ef, 0 denotes averaging over 
all angles and n is the aspect ratio of our Hall bar. 
Measurements are shown after cooling the sample suffi¬ 
ciently below the Neel temperature (to 200 K) and then 
applying a rotating 2 T field. The curves show a har¬ 
monic sm2(pB (cos2(/9b) dependence on the field-angle 
and the amplitudes of AMR\\ and AMR± scale with the 
Hall bar aspect ratio. Note that the crystalline AMR 
contribution is negligible in our case (SI Sec. HI). This 
phenomenology is reminiscent of common AMR traces in 
ferromagnets in applied saturating magnetic fields. 

For comparison, we show in Figs. [^,b also 2 T AMR 
curves measured at a low temperature (5 K). The cor¬ 
responding traces are anharmonic, show history depen¬ 
dence, have smaller magnitudes than at 200 K, and 
AMR\\ and AMR_i do not scale with the Hall bar as¬ 
pect ratio. The MnTe antiferromagnet becomes stiffer at 
5 K and the 2 T field causes only partial reorientation of 
the spin-axes in the domains around their zero-field di¬ 
rection, reminiscent of ferromagnets in weak fields below 
the saturation field. 

In Fig. we show AMR measurements in which, for 
each point, we first heated the sample above the Neel 
temperature (to 350 K) and then field cooled (with Hpc) 
down to 200 K in a 2 T field of a fixed angle (Pb,fc and 
measured the corresponding resistance. To obtain the 
data shown in Fig. Hi we continued with the field-cooling 
down to 5 K then removed the field and took zero-field 
resistance measurements again at 200 K. (Corresponding 
data for other read-out temperatures are shown in SI 
Fig. S6a.) Remarkably, we observe similar AMR traces 
in the two panels only the amplitude of the zero-field 
AMR is about a factor 2 smaller than in Fig. with the 
field on. 
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Note that a 200 K zero-field AMR of a comparable am¬ 
plitude to the one seen in Fig. It is also obtained when 
field-cooling from 350 K down to only 200 K. Crossing the 
Neel temperature in the heat-assisted magneto-recording 
helps the efficiency of the writing process. It is not nec¬ 
essary, however. Setting a stable AMR trace of the form 
of Fig. It, only with a factor 2 smaller amplitude, is 
also possible by field-cooling from 200 K down to low 
temperatures and then performing the zero-field read¬ 
out measurement back at 200 K (see SI Fig. S7). This is 
consistent with the ability, seen in Figs. |^,b, to control 
the antiferromagnetic state by the 2 T field even if below, 
but not too far from the Neel temperature. 

We point out that Fig. not only demonstrates a con¬ 
tinuous harmonic-function AMR in an antiferromagnet 
but it also shows a multiple-stability of states in our 
MnTe memory device. Next we explore how the am¬ 
plitude of the memory read-out signal depends on the 
strength of the writing field and test the limits for eras¬ 
ing the multiple-stable states by turning the magnetic 
field back on. In Fig. we plot the AMR amplitude, 
defined as ^xy(45 ^)-^xy(- 45 obtained when setting 

the states by field-cooling at (fB.FC A 45° from 350 K 
down to 5 K and measuring the resistance during the 
temperature down-sweep with the field on. Fig.[^ shows 
zero-field AMR amplitudes measured during the subse¬ 
quent up-sweep in temperature. 

For the range of writing field amplitudes from 0.5 to 
2 T we were able to set the multiple-stable memory states 
with the zero-field AMR of the harmonic form seen in 
Fig.|2]i. The signal disappears for all traces at ^ 285 K, 
which is consistent with the Neel temperature of thin a- 
MnTe filmP^ and confirms the antiferromagnetic origin 
of the read out signal. The amplitude of the AMR signal 
scales with the amplitude of the writing field, apart from 
the region between 0.5 and 1 T where the AMR signal 
changes sign, and does not saturate at the maximum ap¬ 
plied writing field of 2 T. Before discussing the origin of 
these features we complete in Figs. H ,d the description 
of experimental data, namely of the measurements test¬ 
ing the robustness of the multiple-stable memory states 
under magnetic field perturbations. 

In Fig. li we replot the zero-field AMR amplitude ob¬ 
tained from the temperature up-sweeps after setting the 
states in 2 T writing fields. We compare the trace with 
analogous measurements which only differ in an addi¬ 
tional magnetic field exposure of the memory at 5 K 
after the writing and before starting the temperature 
up-sweep read-out measurement. The additional expo¬ 
sure comprises application of rotating in-plane and out- 
of-plane magnetic fields of a 1 or 2 T amplitude. We see 
that neither of these two amplitudes of the disturbing 
field applied at any angle is sufficient to fully erase the 
memory. After the exposure, the multipe-stable states 
maintain their characteristic harmonic AMR form of the 
read-out signal (see SI Fig. S6) with only the amplitude 
being partially reduced, as shown in Fig. H- 

We performed similar attempts to erase our antiferro¬ 


magnetic memory at 200 K. The results, shown in Fig. p, 
illustrate that at this elevated temperature the 2 T mag¬ 
netic field is sufficient to fully erase the memory states. 
This is consistent with measurements in Figs.[^,b in the 
rotating 2 T field at 200 K, where the observed AMR 
traces had a character of AMR in ferromagnets under 
saturating magnetic fields. AIT field, on the other 
hand, is not sufficient to erase the memory at 200 K; it 
only reduces the amplitude of the read-out signal, as seen 
in Fig. j^. 

In Fig. we illustrate in more detail the stability of 
our antiferromagnetic memory in fields, which are insuf¬ 
ficient to erase it. We explore how the states set by 
cooling in 2 T writing fields of angles (fB,FC = ±45°, 
corresponding to the extrema in the AMR± read-out sig¬ 
nal, are disturbed at 5 K by a 1 T field rotating in the 
sample plane. By taking the resistance measurements 
with the field on we observe a partial reorientation of 
the antiferromagnetic spin-axis, reflected in the varying 
resistance signal. The variations are, however, smaller 
than the difference between the zero-field resistances of 
the two extrema. Moreover, if at each given angle of 
the disturbing magnetic held we remove the held and re¬ 
peat the resistance measurement, we see that the original 
zero-field state almost fully recovers. 

An additional proof of only partial reorientations of 
moments in the antiferromagnetic domains by the ap¬ 
plied rotating held at 5 K is given in Fig. S8 in the SI. 
We observe that the ratio of the amplitudes of AMR\\ 
and AMR± depends on the previous fleld-cooling proto¬ 
col and generally does not scale with the Hall bar aspect 
ratio. This is in contrast to measurements at 200 K in 
the 2 T rotating held (Fig. P ,b) where the AMR shows 
the scaling with the Hall bar aspect ratio, independent 
of the history. 

We now proceed to the discussion of the origin 
of the observed multiple-stable aniferromagnetic mem¬ 
ory states and the corresponding zero-field harmonic- 
function AMR. When cooling the system in an applied 
magnetic held, the domain distribution freezes at temper¬ 
atures sufficiently below the Neel temperature, as illus¬ 
trated in Fig. P- The freezing occurs when the thermal 
fluctuations no longer allow the different domains to flip 
their spin-axis from one to the other metastable direc¬ 
tion. 

The writing magnetic held applied at a certain angle 
during the cooling makes the easy axis closest to the held 
normal more favorable (see sketches in Fig. p) . Due 
to the finite temperature, however, the other two easy 
axis directions can be also populated with the statistics 
given by the Boltzmann distribution and by the relative 
strength of the terms contributing to the total energy 
of the system for the given domain size, namely the ex¬ 
change energy, the anisotropy energy, and the Zeeman 
energy. 

When the Zeeman term dominates, i. e. for very large 
magnetic fields or domain sizes, one of the three easy 
directions would be strongly favored within a 60° in- 
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terval of the writing field angles, and the other direc¬ 
tions for the other respective 60° intervals. The read¬ 
out AMR signal in this case would take the form of a 
step-like trace with three distinct memory states. Upon 
reducing the Zeeman contribution, the AMR trace under¬ 
goes a transition into a continuous function whose am¬ 
plitude gradually decreases and whose shape approaches 
the harmonic b,fc function of the writing field an¬ 

gle. Since our measured signals plotted in Fig. show a 
weak but clearly detectable deviation from the harmonic 
sin 2(pB^FC form we can infer the typical domain size from 
fitting the model to the measured data. 

The distribution of the three different easy-axis do¬ 
mains and the corresponding AMR signals are calculated 
using the antiferromagnetic Stoner-Wohlfarth modeP^. 
The results of the modelling, plotted in Fig. [^, repro¬ 
duce well the experimental behavior. In the model, the 
exchange energy was estimated from the Neel tempera¬ 
ture and the mag netic moment found in neutron diffrac¬ 
tion experimentThe anisotropy energy was then 
inferred assuming a spin-flop field of ^ 2 T. (For more 
details on the model see SI Sec. V.) 

From the fitting we obtain that all three domains are 
at least partially occupied with the highest inequality of 
2 : 1:1 and 2 : 2:1 for the cases when the writing field is 
perpendicular or parallel to one of the three easy axes, 
respectively. The moments easy orientation is found to 
be along (lOlO) directions as indicated in Fig. [^. The 
domain size we find in the fitting is approximately 20 nm, 
which is in a good agreement with the mosaic block size 
obtained from the X-ray diffraction results. From this we 
conclude that the magnetic domain size is limited by the 
structural block size and does not change significantly 
with the writing conditions. 

Under the assumption of a constant domain size we use 
the model to also estimate the dependence of the ampli¬ 
tude of the zero field read-out signal on the strength of 
the writing field. As shown in Fig.[^, the model predicts 
an initial quadratic increase of the read-out AMR sig¬ 
nal followed by a tendency towards saturation at higher 
fields. It broadly agrees with the measured data apart 
from the negative AMR signal seen in the experiment at 
low writing fields. We attribute this negative signal to 
domain walls which may have some uncompensated mo¬ 
ments. Unlike the compensated moments in the antifer¬ 
romagnetic domains which tend to align perpendicular to 
the applied field, these uncompensated moments prefer 
a parallel alignment. Since the uncompensated moments 
only lead to a linear gain in the Zeeman energy the ef¬ 
fect is diminished at larger fields by the canting of the 
moments in the antiferromagnetic domains which yields 
an energy gain quadratic in the applied magnetic field. 

To summarize, we have reported spintronic memory 
functionalities in an antiferromagnetic counterpart of 
common II-VI compound semiconductors. Favorable 
magnetic characteristics of our a-MnTe epilayers allowed 
us, on one hand, to evidence the direct analogy between 
AMR in antiferromagnets and ferromagnets, both in its 


basic phenomenology and in the utility as an electrical 
detection tool of the ordered magnetic moments. On the 
other hand, our work also highlights the unique poten¬ 
tial of antiferromagnets for spintronics. In our MnTe de¬ 
vices we demonstrated that the harmonic-function AMR 
persists even after removing the setting magnetic field 
and that the multiple-stable memory states cannot be 
erased by strong magnetic field perturbations when suf¬ 
ficiently below the antiferromagnetic transition temper¬ 
ature. These magnetic memory characteristics are un¬ 
precedented in ferromagnets. 


METHODS 

Sample preparation. o-MnTe was grown by molec¬ 
ular beam epitaxy on (111) oriented InP and SrF 2 sub¬ 
strates using elemental Mn and Te sources and substrates 
temperatures in the range of 370 to 450° C. 2 D growth 
was achieved in both cases as judged from the streaked 
RHEED patterns observed during growth. The orienta¬ 
tion of our layers on the substrate is (0001)[10l0]MnTe|| 
(lll)[112]sub- For structural characterization we carried 
out cross-sectional scanning transmission electron mi¬ 
croscopy (STEM) in the high angle annular dark field 
(HAADE) imaging mode and x-ray diffraction studies 
with CuKo^i radiation. STEM-HAADE images were 
acquired with a NION UltraSTEM, equipped with a 
NION aberration corrector and operated at 100 kV. In 
the images the intensity scales approximately with the 
atomic number (Z) squared.!^ Cross sectional STEM 
specimens were prepared using a EEI Nova 200 Dual- 
Beam SEM/EIB focused ion beam. 

Magnetometry measurements were performed on 
2000 nm thick MnTe films with area of ^ 11 mm^ 
and 500 /im thick InP substrate in a Quantum Design 
SQUID magnetometer using reciprocating sample option 
(RSO) for increased measurement sensitivity. Magnetic 
field sweeps were recorded for various crystallographic 
directions while sweeping from negative to positive and 
vice versa showing no hysterisis. Erom all magnetic field 
sweeps a negative slope obtained from the fit to the hard- 
axis out-of-plane sweep seen in the inset in Eig. was 
subtracted to correct for signal coming primarily from 
the diamagnetic substrate. Temperature dependent sus¬ 
ceptibility measurements were taken in magnetic field of 
0.5 T. 

Optics The absorption coefficient and thus band 
gap properties were determined by fitting the transmis¬ 
sion spectra (measured at 300 K) with an analytical 
expressiorP^ for the transmission of a thin absorbing film 
on a transparent substrate. 

Processing and transport measurements Eor 
transport measurements Hall bars were defined in 50 nm 
thick films grown on InP along the [lOlO] direction 
by electron-beam lithography and wet etching using a 
H 3 PO 4 /H 2 O 2 /H 2 O mixture as etchant. Contacts were 
made by soldering with pure In and 4-point measure- 
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ment technique was employed to avoid possible contribu¬ 
tion of contact resistances. Complementary Corbino-disk 
devices were fabricated by deposition of An rings. The 
transport measurements were performed in an Oxford 
Instruments vector field cryostat capable of 2 T in arbi¬ 
trary directions and up to 6 T in the current direction. A 
constant DC current density of 2000 A/cm^ was applied 
during all measurements. 
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FIG. 1. MnTe thin films properties, a) Cross sectional high resolution high angle annular dark field (HAADF) image 
taken along the [1210] zone axis resolving the individual atomic columns of Mn (magenta) and Te (green). The inset shows the 
atomic positions including the Mn magnetic moments in the ordered antiferromagnetic state and the unit cell size is indicated 
by a white line, b) Absorption coefficient of MnTe in the mid-infrared and visible spectral range extracted from transmission 
measurements for three different film thicknesses, c) SQUID magnetic field sweeps for different field directions at 5 K. The 
slope of the out-of-plane measurement was subtracted from all the measured traces. An inset shows the as measured moments 
dominated by the diamagnetism of the substrate, d) Temperature dependent resistance measurement showing a peak around 
the transition temperature. An inset shows the Hall bar geometry and the definition of the magnetic field angle (fB- 
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FIG. 2. Antiferromagnetic anisotropic magnetoresistance and multiple-stable memory. a,b) Transverse (red) and 
longitudinal (blue) AMR measurements at 200 K and rotating in-plane 2 T field. Analogous measurements at 5 K are shown 
in grey. Arrow in panel a) indicates the initial angle and the direction of rotation, c) Transverse AMR measured at 200 K 
after cooling from 350 K in a 2 T magnetic field applied at an angle ipB,FC and with the field kept on. d) Zero-field transverse 
AMR obtained after continuing with the field-cooling down to 5 K then removing the field and taking zero-field resistance 
measurements at 200 K. The red line shows the fit of the model calculations described in the text, while the grey line is a 
~ sin 2(pB,FC fit. 
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FIG. 3. Dependence on writing field and robustness against field perturbations, a) Temperature and writing field- 
strength dependence of the AMR amplitude for the temperature down-sweep with the writing field kept on. b) Temperature and 
writing field-strength dependence of the zero-field AMR amplitude for the temperature up-sweep, c) Zero-field AMR amplitudes 
for the temperature up-sweep after exposing the memory at 5 K to 1 T (green) and 2 T (red) perturbing fields rotated both 
in-plane and out-of-plane. Blue curve represents the unperturbed reference measurement obtained with Bpc = 2 T. d) Same as 
c) for the field perturbations applied at 200 K. e) The stability of two memory states set by the heat-assisted magneto-recording 
from 350 K with a 2 T writing field applied at angles (fB.FC — ±45° tested at 5 K by a rotating 1 T field. At every angle (fB 
of the perturbing 1 T field, red/green curves correspond to the read-out resistance measurements with the field on while black 
lines are obtain after removing the perturbing field at each if b ■ 
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FIG. 4. Modeling of the multiple-stable states, a) Schematics of setting the multiple-stable states, each with a 
distinct domain distribution, by the heat-assisted magneto-recording across the Neel temperature, b) Schematics of the two 
antiferromagnetic spin sublattices which tend to align perpendicular to the applied magnetic field and to cant towards the field; 
schematics of the three in-plane easy axes in the hexagonal o-MnTe. c) Modelling of the distribution of the three domains 
in the multiple-stable memory states and of the corresponding transverse AMR signal for different size of the domains, d) 
Comparison of the model prediction (black line) and experimental dependence at 200 K (red points) of the read-out AMR 
amplitude on the strength of the writing magnetic held. 
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Supplementary Informations 

Multiple-stable anisotropic magnetoresistance memory in an¬ 
tiferromagnetic MnTe 


I. MNTE THIN FILMS - GROWTH AND BASIC PROPERTIES 


MnTe thin films were grown on InP(lll) and SrF 2 (lll) substrates for transport and optical investigations, respec¬ 
tively. Growth was performed by molecular beam epitaxy using elemental sources at a substrate temperature of 370 
to 450° C. Using X-ray diffraction we find that the thin films grow in the hexagonal NiAs bulk phase (a-MuTe) and 
that no other phases are present (Fig. [sT|i). In magnetometry measurements, although the signal is dominated by 
the temperature independent diamagnetism of the substrate, we are able to detect the transition temperature around 
^ 300 K (Fig.|Sl|3) from a peak in the susceptibility. Furthermore we find that the magnetic moments are oriented in 
the sample plane as evident from the drop in susceptibility which is more pronounced for an in-plane magnetic field, 
equivalent to the report in Ref. 21 (main text). Temperature dependent diffraction experiments reveal an anomaly of 
the out of plane lattice parameter at the transition temperatur^^ (Fig. S2), confirming the transition temperature 
found by the magnetometry measurements. 


The semiconducting properties of our a-MnTe thin films are clearly seen in the Hall measurement shown in Fig. SI 


which shows a linear dependence of the Hall voltage on an out of plane applied magnetic field. We find a hole carrier 


density of p = 6 x 10 


18 3 


with a mobility of /i = 43 cm^/Vs at T = 5 K. 
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FIG. SI. MnTe thin film properties, a) Xray diffraction scan of a 50 nm thick MnTe film grown on InP(lll) showing 
only phases of the hexagonal bulk o-MnTe phase which is sketched in the inset, b) Magnetic susceptibility of o-MnTe thin 
films measured for in-plane and out of plane field orientation. The drop of the susceptibility in the in-plane configuration is a 
result of the in-plane moment orientation, c) Hall measurement in order to determine the carrier density and mobility. 
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FIG. S2. temperature dependent c lattice parameter, a) The out of plane lattice parameter shows a anomaly in its 
thermal expansion properties at the AFM transition temperature, b) In the derivative of the lattice parameter the anomaly is 
even more pronounced. Measurements were performed on a 2000 nm thick a-MnTe grown on InP. The Neel temperature Tn 
is indicated by arrows. 


A. Epitaxial orientation determined by X-ray diffraction 

Even though the mismatch between a-MnTe and the InP (111) or SrF 2 (lll) substrates is below 1% the films grow 
relaxed for the investigated thickness range of 50 to 2000 nm. From the peak positions in the reciprocal space maps 
we find lattice parameters of a = 4.1708 A c = 6.6860 A for 50 nm MnTe grown on InP. Further we are able to 
determine the in-plane epitaxial orientation of the hexagonal a-MuTe lattice on the cubic substrate. For both types 
of substrates the orientation of the c-planes (0001)1^ is parallel to the (111) planes of the cubic substrate and the 
in-plane [lOlO] direction of a-MnTe corresponds to the [112] direction of the substrate. Figure shows reciprocal 
space maps measured for a-MnTe thin films grown on both types of substrate. With a mosaic block modeP^ (dashed 
contour lines in Fig. |S3[ ) we are able to describe the peak shape of the diffraction signal for both types of samples 
(grown on InP(lll) or SrF 2 (lll)). For the sample used in the transport measurements (grown on InP) the lateral 
and vertical blocksize is found equal around 25 ± 5 nm and a Gaussian rotation distribution with standard deviation of 
0.25 ±0.1 degree was used for the simulations. The sample used in optics (grown on SrF 2 ) can be described by blocks 
of same lateral size, however, with 50 ± 2 nm vertical block size and Gaussian rotation distribution with standard 
deviation of 0.15 ± 0.05 degree. 
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FIG. S3. X-ray diffraction reciprocal space map measurements. The color plot shows the experimental data with 
labeled peaks from the SrF 2 /InP substrate and o-MnTe film. The shown measurements are for 200 nm MnTe on SrF 2 in panels 
(a-c) and 50 nm on InP in (d-f). Dashed contour lines mark the fit with the mosaic block model. 


II. OPTICAL CHARACTERIZATION 

We determine the spectral dependency of the absorption coefficient for samples with different film thickness and 
find low absorption in a wide spectral range as typical for a semiconductor. Below a 700-800 nm all the films show 
an onset of strong absorption due to excitation of carriers across the band gap. The used samples were grown on 
SrF 2 , which is transparent in the full visible spectral range and up to 10 /im in the infrared (Fig. [S4^b) ruling out 
any influence of the substrate. Due to the film thickness we observe Fabry-Perot oscillations in the transmission data 
at wavelength above the onset of the strong absorption. By modeling (see Ref. 26 in main text) these oscillations 
considering the sample thickness we extract the thickness independent absorption coefficient from the measurement. 
The spectral dependence in Fig. lb (main text) is typical for an indirect bandgap material. The size of the band gap 
at room-temperature is found as Eg = 1.46 ± 0.10 eV. 
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FIG. S4. Optical characterization of MnTe by transmission measurements. Transmission spectra of the o-MnTe 
thin films grown on SrF 2 (lll) in the visible (left) and mid-infrared (right) spectral regions. Data for samples with various 
thicknesses allow us to extract a value of absorption coefficient shown in the main text. 


III. MAGNETIC FIELD THRESHOLD AT ELEVATED TEMPERATURES 

From the magnetometry data shown in the main text we conclude that the low temperature onset of a sizable mag¬ 
netic moment occurs at a ^ 2 T in-plane field which may be associated with the spin-flop field in the antiferromagnetic 
domains. For higher temperatures the magnetometry measurements are shown in Fig.|S5^. These measurements were 
also correlated with in-plane field sweeps in a Corbino geometry. The longitudinal magnetoresistance shows a similar 
change of slope at low temperatures. This feature remains clearly visible at higher temperatures and allows us to 
conclude that the critical field needed for a reorientation of the moments decreases at higher temperatures. At 200 K 
the kink is sufficiently below 2 T to allow manipulation of the Neel vector by application of a 2 T field. From the 
Corbino geometry we further can estimate the size of the crystalline anisotropic magneto-resistance (AMR), which 
was found to be below 0.03% at 200 K and 2 T. It is therefore a negligible contribution in the transport data presented 
in the main text. 
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FIG. S5. Magnetometry measurements and correlation with resistance at elevated temperatures, a) Magnetic 
moment in in-plane magnetic field sweeps recorded at different temperatures. The temperature independent diamagnetic 
contribution of the substrate was subtracted as explained in the main text. The kink around 2 T gradually disappears at higher 
temperatures, b) Transport measurements in Gorbino geometry at various temperatures show that an analogous kink can be 
observed in the resistance which via the crystalline AMR effect is sensitive to the orientation of the magnetic moments. 
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IV. STABILITY OF THE FROZEN STATES 


In the main text we show the zero field AMR signal at 200 K (Fig. 2d) as it results after field cooling from above Neel 
temperature to 5 K. Figure [S 6 ^ shows the equivalent data for different temperatures. The AMR signal is maximum at 
150 K and reduces to zero at the magnetic order transition. To test the stability of this signal with respect to magnetic 
field at low temperature we performed field rotations at 5 K. Although the zero field AMR signal is decreased as can 
be seen in Fig. [ 86 ) 3 , c it can not be fully erased by neither 1 T nor 2 T rotations. This is in contrast to equivalent 
rotations at 200 K, where a rotation of a 2 T field (above the spin-flop field) is able to destroy the frozen state as seen 
in main text Fig. 3d. 

The fact that with 2 T at 200 K we are in a different regime as at 5 K is also seen in the ratio of the AMR\\ to 
AMRj_ amplitudes. Only at 200 K and 2 T the respective amplitudes scale with the Hall bar aspect ratio. In the 
other cases (5 K: 1 and 2 T, and 200 KIT) the data show a dependence on the field cooling direction and strongly 
deviate from the ideal ratio of 1 . 



(pB,Fc (deg) 

FIG. S6. Stability of the AMR signal against field rotations, a) The zero field AMR signal induced by field cooling 
in various directions is shown for different temperatures. Below 300 K there is clear semi harmonic signal due to the frozen 
moment orientation. The same signal with equal symmetry is found in a measurement after the magnetic field of a) 1 T and 
b) 2 T magnitude was rotated in three perpendicular planes at 5 K. Although the amplitude of the AMR is slightly reduced 
these field rotations can not destroy the frozen state. 
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FIG. S7. AMR|| to AMR^ ratio for different field cooling directions. The ratio of the AMR signal extracted from 
the longitudinal and perpendicular contacts shows strong dependence on the field cooling direction at 5 K. At 200 K and 2 T 
the ratio is independent of the history of the sample indicating that 2 T are above the critical field for reorientation of the 
moments. Dashed circles mark the data points extracted from the curves shown in Fig. 2a,b (main text). 


Due to the fact that at 200 K we are able to overcome the spin-flop field with a 2 T in-plane field we are also 
able to perform magneto-recording when cooling in field from 200 to 5 K. Corresponding data are shown in Fig. [S8| 
and clearly demonstrate that at 200 K AMR traces are found in zero field although the Neel temperature was never 
overcome. 
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FIG. S8. Magneto-recording by cooling from 200 to 5 K. AMR± data recorded in zero field at 200 K after field cooling 
from 200 to 5 K in a 2 T applied field. 
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V. MULTIPLE-STABILITY AND ANISTROPIC MAGNETORESISTENCE MODEL 

In order to describe the domain population after the field cooling and the resulting AMR signal we use the Stoner- 
Wohlfarth model. We apply it to each magnetic sub-lattice of a-MuTe in order to describe the anisotropy of the 
material and we introduce a term accounting for the exchange energy. We restrict ourself to a 2D description which 
is sufficient to describe the field cooling in different in-plane magnetic fields. Every magnetic domain consists of two 
equal ferromagnetic sub-lattices which have opposing magnetization vectors Mi^ 2 - When a magnetic field B is applied 
in the plane of the sample the magnetic moments might cant with respect to each other and thereby tilt away from 
the easy axis direction. The relevant angles are defined with respect to the current or easy axis direction as sketched 
in Fig. 

The total energy E per sample volume U of a single domain is given by the sum of exchange, Zeeman, and 
magneto-crystalline anisotropy energies: 


E/V = JexMi • M2 — B • (Ml + M2) + Emae(Mi) + Emae(M2) 


( 1 ) 


where Jex is the exchange constant, Emae is the magneto-crystalline anisotropy funtion, and a hat (” "”) denotes a 
unit vector. Using the angles denoted in Fig. S9 and the exchange field defined as = Je^jM with M = |Mi^ 2 | 
this can be rewritten as 


E/V = - cos('0i -h '02) - mb [cos {ipB - 2 pi)) - cos {^pB - {a- -02))] 

+ £’mAe('0i) + ^MAe('02) (2) 

The exchange field we obtain from the Neel temperature as Bq^ = where we use the Boltzmann constant 

kB and the Bohr magnetron /tb- The magneto-crystalline anisotropy energy density we express as 

^MAE (0) = ^MAE sin^ (30), (3) 

in order to describe the six-fold in-plane symmetry of the hexagonal material. Amae has units of an energy density 
and can be obtained from the spin flop field Bsf and the exchange field by Amae = 12M5|p/5ex- 

From the descibed model we can obtain the energy gain AEi for the three different easy axis directions ai when 
a magnetic field is applied at a certain angle. For this we find the arrangements of the moments by minimizing the 
energy for every field direction. Using the obtained energy gain for a certain sample volume V in comparison with 
the thermal energy at the moment of the freezing in the Boltzmann distribution we obtain the relative population wi 
of the different domains. 


pAEi/(/cBTF) 

Wi = —o- (4) 

y-3 AE,/(/cbTf) ^ ^ 

This population numbers we use to calculate the transversal and longitudinal AMR signal which is proportional to 
Y/^^iWiSm{2ai) and cos(2(a0, respectively. Calculating the expected zero field AMR for different domain 

sizes one can see how the field cooling efficiency is increasing for bigger domains leading to stronger AMR (Fig. 4c 
in the main text). When the domains are getting big enough so that almost only one domain is populated the AMR 
signal is getting step-like. The results of our calculations in comparison with the experimental data, which show a 
transitional behavior between the harmonic AMR and the step-like signal, are shown in the main text Fig. 2d. By 


B 



FIG. S9. Sketch of the parameters in the Stoner—Wohlfarth model. The angles of the magnetic moments of the two 
sub-lattices deviate slightly from the easy axis direction upon the application of an magnetic field. For the modeling relevant 
vectors and angles are indicated. 
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adjusting the domain size and easy axis direction as well as amplitude of the AMR a nearly perfect agreement with 
the experimental shape could be found. The easy axis direction is found to be along the (lOlO) directions and the 
domain size was found around 20 nm similar to the domain size determined by the mosaic block model. The variation 
of the AMR amplitude with field cooling strength also qualitatively agrees with the experimental observations. 
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